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ABSTRACT 


The  effort  here  was  concerned  with  initiating  a  long 
term  experimental  study  to  develop  the  method  of  laser 
light  scattering  as  a  plasma  diagnostic  technique.   This 
was  done  with  a  ruby  laser  and  a  grating  spectometer  cou- 
pled with  sensitive  detectors,  at  the  steady  state  plasma 
facility  at  the  Naval  Postgraduate  School.   Although  at- 
tempts to  detect  the  extremely  low  levels  of  scattered 
light  failed,  excellent  progress  was  made  in  determining 
existing  apparatus  limitations  and  in  recognizing  which 
system  parameter  may  be  adjusted  in  order  to  detect  the 
scattered  light;   this  will  enable  the  future  experimenter 
to  progress  with  a  minimum  of  overlap. 
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I.   INTRODUCTION 

This  thesis  is  intended  as  a  detailed  report  of  the 
beginning  of  a  long-term  experimental  effort  to  develop 
Thomson  scattering  of  laser  beams  as  a  plasma  diagnostic 
technique.   The  technique  employs  the  particle-like  behav- 
ior of  photons  in  scattering  off  plasma  electrons.   In 
general  the  background  radiation  from  an  energetic  plasma 
dominates  making  it  impossible  to  distinguish  the  scatter- 
ed light.   This  problem  may  be  somewhat  alleviated  by  the 
use  of  an  extremely  intense,  sharply  tuned  source  in  con- 
junction with  narrow  band  detectors.   The  high  power  ruby 
laser  is  such  a  source  and  a  grating  spec trograph -photo 
multiplier  tube  arrangement  is  such  a  detector.   Even  with 
these  improvements,  a  manifold  of  problems  and  complica- 
tions yet  exist,  not  the  least  of  which  is  careful  planning 
of  the  geometry  of  the  experiments  and  associated  alignment 
of  apparatus.   Previous  experiments  have  shown  that  laser 

scattering  from  electron  beams  of  energy  approximately  2 

9    -3 
kev  and  density  approximately  5  x  10   cm   yielded  scatter- 

— 18 
ing  intensities  of  only  about  10    of  the  laser  output. 

Though  electron  densities  in  the  plasma  under  study  are 
significantly  higher,  this  gives  some  idea  as  to  the  prob- 
lems presented  in  discerning  the  scattered  light. 


It  will  be  shown  that  the  spectral  distribution  of  the 
scattered  light  is  in  part  determined  by  a  dimensionless 
coherence  parameter  proportional  to  the  ratio  of  incident 
light  wavelength  to  the  Debye  length; 

o 

when  Oi   is  small,  that  is,  when  the  wavelength  of  incident 
light  is  small  compared  with  the  Debye  length  the  cross 
section  is  the  appropriate  Thomson  cross  section  and  the 
shape  of  the  scattered  spectrum  is  a  Gaussian;   the  width 
is  the  Doppler  width  associated  with  electron  motion  and 
as  such  is  a  measure  of  electron  temperature. 

When  a,   is  large  (G?  »  1)  the  spectrum  has  a  central 
line  of  the  incident  light  frequency,  L£     ,  between  two  sat- 
ellite peaks  at  frequencies  0?  +  to  where  a?  is  the  plasma 
frequency.   In  this  case  the  Doppler  broadening  of  the  cen- 
tral peak  is  approximately  proportional  to  the  ion  motion 
and  gives  information  concerning  ion  temperatures. 

The  research  begun  in  the  plasma  laboratory  of  the 
Naval  Postgraduate  School  has  been  under  the  cognizance  of 
Professor  F.  Schwirzke  and  is  financed  under  Naval  Air  Sys- 
tem Command  Research  Grant.   The  accomplishments  made  thus 
far  and  described  herein  will  enable  to  the  future  experi- 
menter to  progress  with  a  minimum  of  overlap.   The  major 
obstacle  yet  to  be  overcome,  and  indeed  the  main  problem 
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in  the  development  of  this  diagnostic  method  as  a  whole,  is 
detecting  the  scattered  light  and  distinguishing  it  from 
natural  background  radiation.   Once  this  is  done,  the  tech- 
nique may  be  rapidly  and  effectively  developed. 

The  diagnostic  technique  once  developed,  it  may  be 
employed  to  investigate  turbulence  in  a  laser  produced 
plasma. 
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II.   METHOD  AND  THEORY 

A.    ELECTROMAGNETIC  THEORY  OF  LIGHT  SCATTERING 

Following  directly  from  the  equation  of  motion  the  am- 
plitude of  an  electron  oscillation  caused  by  the  electric 
field  E_'  of  the  incident  electromagnetic  wave  of  frequency 

Ui     is; 
o 

s(t,  -  slm  =  4^  (1) 

m  CO       4ff  mc 
o 

Multiplying  the  displacement  by  the  electron  charge  e 
the  induced  dipole  moment  is  obtained; 


-e  X  E        , 
p  =  -  es_   =   — - — y-   =  cE  (2) 

4ff  mc 


For  the  laser  light  scattering  to  be  dealt  with  here, 
is  a  plane  monochromatic  wave  cut  off  exponentially  at 
time  T  =  1/y  ,  the  pulsetime  constant. 

That  is, 

E'(r,t)  =  E'  expCiCk'-r-f  (050+iy)t]}  (3) 

p(r,t)   =  eEgexp{i[k'  -r  +  (0>o  +  iy)  t]  }  (4) 

It  is  this  time  dependent  Dipole  moment  (4)  which  produces 
the  scattered  radiation. 


10 


The  electric  field  E."  (r,t)  scattered  by  the  j    electron 
at  position  r.  detected  at  R,  from  the  Hertz  formula  for 
dipole  radiation,  is 

E."(f*,t)  =  ^(fr^t)  exp[i(k'.r..  +  k'R.J  (5) 

Rj  =  «j  "  £j  (6) 

where  R.  is  the  distance  between  scattering  center  and  de- 
3 

tec  tor 


f.(R.,t)  is  the  amplitude  of  the  scattered  field  and  is 

-D    ~3  * 

2 

f.(R.,t)=   -^fr    [R.    X     (R.    X    CE'jlexpCiO)  -y)  t]  (7) 


2 2  ) 

477  mc 


£*<&*>*)    =   "   (-%)     R*Cii    x(a.x  E')]exp(iw  t-yt)       (8) 

J       J  mc  _j      J  J 

R. 

where   a . =     zr  ■         is    the   unit   vector   in    the   R.    direction 

:  i  Rj  |  3 

r.  is  measured  from  a  reference  point,  0,  in  the  scattering 
volume.   If  the  dimensions  of  the  scattering  volume  are 
small  compared  with  R.,  then  R.  may  be  replaced  by  ft.,  and 
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f_.  =  f   for  all  electrons  in  the  plasma  volume.   When  deal- 
ing with  the  phase,  however,  more  care  must  be  taken  in 
approximations.   The  following  substitution  may  be  made  for 

r  .  «  ft  . 

-3  ~ 


R.  =  7(ft-r.)2'  =  y(R2_2ft.r.+r>2'  M  ft   -(ft'r^) 

-3         — 3  — d  :     R 


+ 


(9) 


Writing  the  scattered  wave  vector  k"  as  having  approximately 
the  magnitude  of  the  incident  wave  vector  k*  and  direction 
ft,  that  is, 


-  "  k'  I  ft! 


(10) 


(5)  may  be  rewritten  as 


e"  (_ft,  t)  =  f  (ft,  t,X9)exp(ik'ft  +  ik.r.) 


where  k  is  the  scattering  vector.   That  is 


(11) 


4ff  "\  .  Q        ,  .    ,  „    2ft 


"         1  I      4"  \  O  11 

k  =  k  -k  ,   k  =  ^  -y-  )sin  y  ,  k '  «  k   = 


(12) 


The  scattering  angle  6  is  between  k   and  k  .   The  direction 
is  of  course  determined  by  the  incident  laser  beam  and  by 
the  position  of  the  detector. 

In  90°  scattering, 

T 

1 
I 

k, 


\±-#-~* 


da 
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Furthermore  X'?«  X" =  X.   The  wavelength  difference  between 
X'  and  X"  is  usually  small  and  may  be  neglected,  such  that 
the  symmetric  of  the  scattered  spectrum  is  undisturbed;   if 
o>  is  the  frequency  of  the  incident  light  and  o>  is  the  fre- 
quency shift  of  the  scattered  light,  I(o>  +0))«  1(05  -03)  if  no 
drift  velocities  exist. 

The  phase  of  the  scattered  light  depends  on  the  posi- 

ik  •  r  . 
tion  of  the  scattering  center  by  the  factor  e    -1  .   This 

becomes  important  when  the  fields  scattered  by  a  number  of 

centers  are  superimposed  at  _R.   Summing,  over  all  electrons 

in  the  scattering  center,  one  obtains 


N  N  j,    ■_ 

E" (R,  t)  =  ^E.(R, t)  =  f(R, t)e"lk'R  Je  ~      j       (13) 
j=l  j=l 


The  scattered  intensity  is  equal  to  the  time  averaged 

i   "I2  C 
Poynting  vector  |  E    -£=   ;   Introducing  the  incident  light 

intensity  as 

i   i  2 

I    =   I  s  I — £.  (14) 

one  obtains 


c-  i(k-r_.)  2 
1(6, X)  =  i(k)  =  a  (6)    <  \l   e      :  I   >  (15) 


j=l 
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Where  <J(0),  Thomson's  scattering  factor,  is  equal  in  essence 
to  |  f  I;   Introducing  e'  and  e",  the  polarization  vectors 
of  the  incident  and  scattered  light, 

E'  =  e  *|  E '  |  ,   E"  =  e"|E"|   ,   fo-e"  =  0  (16) 


,2 

<T(8)  =  lQs.  =   x 


,   , .  4ff2 
e  •  (ce')  —T- 

X  % 


(17) 


for  linearly  polarized  light;  and  for  natural  incidents 
light, 


mc  2  5* 


where  the  1  in  parentheses  represents  the  contribution 

2 

from  natural  incident  light  and  cos  0  represents  the  con- 
tributed polarized  parallel  to  the  plane  of  scattering.  It 
is  easily  seen  how  (17)  complies  with  this. 

B.    SPECTRAL  INTENSITY  DISTRIBUTION  OF  SCATTERED  LIGHT 

For  the  most  part,  it  is  the  electron  density  fluctu- 
ations in  the  plasma  which  are  responsible  for  the  scatter- 
ing of  light;   that  is,  according  to  Huygen ' s  principle  no 
scattering  can  occur  in  a  strictly  homogeneous  medium.   The 
electron  density  n(r_,  t)  is  a  function  of  position  and  time, 
and  may  be  considered  as  a  continuously  varying  quantity  or 
as  a  superposition  of  delta  functions  as  circumstances 
warrant. 
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Imagine  the  plasma  confined  to  a  large  cubical  volume 

3 
V  =  L   and  write  the  electron  density  as  a  spacial  fourier 

1 
series . 


n(r,t)  =  i  £  n^tje1^  (19) 

q  - 


Where  q  is  the  wave  vector  satisfying  the  cyclic  boundary 
conditions 

qi  =  ^l^  h    =  °'±:L'±2    >  hi  =  x>y>z  (2°) 

The  fourier  coefficient  n  may  be  considered  as  the  ampli- 
tude of  the  longitudinal  plasma  wave  with  wave  vectors  q 
and  may  be  written 


n   =   n(r,t)  exp  (-iq.rjdr  (21) 

-    V 

for  the  microscopic  representation  and 

N 

nq  =  L   exp(-iq-rj)  (22) 

j=l 

for  the  microscopic  representation. 

Usually,  longitudinal  oscillations  in  plasmas  have  com- 
plicated time  dependencies  and  dispersion  relations;  further- 
more a  number  of  frequencies  may  be  present.   The  probability 


1 
As  done  in  Light  Scattering  and  Drift  Instability  in 
a  Plasma  by  Otto  Theimer.  ltc. 
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of  a  certain  frequency  occurring  is  related  to  the  mean 
square  of  the  amplitude  n  (to)  of  a  wave  with  wave  vector 
q  and  frequency  05  .   A  long  lived,  relatively  undamped 
oscillation  will  have  a  sharp  maximum  in  n  (tf)  .   A  con- 
tinuous function  of  the  frequency  n  (to)  is  the  Fourier 
transform  of  n_(t)  with  respect  to  time  and  may  be  written 
in  integral  form; 


nq(05)  =  J  n  (t)  exp(-icot)dt  (23) 


But  the  actual  density  fluctuation  is  not  by  itself  of  pri- 
mary concern  in  scattering  theory;   rather,  the  amplitude 
of  the  polarization  waves  P  ( t)  induced  by  the  incident  ra- 
diation is  of  interest.  [Not  to  be  confused  with  the  dipole 
moments  of  equation  (2) ] 

P  (t)  =  €    n  (t)E'  (t)  (24) 

Where  c  is  given  by  (7a)  and  is  the  polarization  of  the 
free  electrons.   Again,  E  (t)  is  the  electric  field  of  the 
incident  laser  pulse  and  has  the  time  dependence 

E'  (t)  =  E1  expC  (io)  -y)  t]  ,    t  >  0 

(25) 
E  (t)  =  0  ,    t  <  0 

That  is,  the  laser  is  fired  at  t  =  0  and  is  essentially 
overafter  T  =  —  .   Thus  it  is  seen  that  the  polarization 

y 

waves  have  a  frequency  dependence  (to  +co)  the  amplitude  of 
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which  decreases  exponentially  with  the  laser  pulse,  as 
exp(-yt).   That  scattered  light  will  of  course  have  the 
same  time  dependence,  and  the  scattering  waves  may  be  iden- 
tified as  a  frequency  shift  in  the  scattered  spectrum. 
Thus  the  finite  laser  pulse  may  be  attributed  (at 

least  mathematically)  to  the  plasma  waves  and  n  (to)  may  be 

~i 

written  as 


n  (a))  = 


n  (t)  exp  [-(iw+  y)  tldt 
o  -i 


(26) 


Considering  the  total  scattered  intensity,  the  time  aver- 
age value  of  the  function  n,  (t)  is  taken 


<h,<t>|2>  = 


nk(t)|   exp(-2yt)dt 


exp(-2yt)dt 


o 


00 

=  2y    [^(t)  |   exp(-2yt)dt 


(27) 


By  Parseval ■ s  Theorem     from  (26)  and  (27)  it  follows 
that 


CO 

<k<fc)l2>  - ?  J     k(a»  I2  da 


(28) 


_00       _ 


Clearly,  additively  combining  the  square  of  the  amplitudes 
(^•)  I  n,  (W)|    gives  the  square  of  the  amplitude  ^|  n,  (t)|  /, 


In  Bronstein  and  La.  Semgndjajew,  Taschenbuch  d.  Mate- 
matik,Verlag  Harri  Deutsch, Frankfurt/M. , 1965, va5  pp. 404-474 
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Thus    each  wave   of  wave   vector  k   and   frequency   CC  contributes 
to    the    total    scattered   intensity    by   an   amount   given  by 

I(W  ,k)dO>  =   (7(6)    (|  )  1^(05)  I2da)  [  ?«S|£]  (29) 

And  as  already  stated,  each  frequency  shows  up  as  a  fre- 
quency shift  in  the  scattered  spectrum. 

So  it  has  been  established  that  each  electron  becomes 
Hertz  oscillator,  when  induced  to  oscillate  by  the  inci- 
dent light,  and  itself  emits  the  light  which  is  the  scatter- 
ed radiation.   If  all  electrons  oscillate  in  phase  with  the 
incident  light,  the  emitted  (scattered)  radiation  will  be 
coherent  in  time.   If,  however,  the  electrons  have  complete- 
ly arbitrary  spacial  distribution,  coherence  is  destroyed. 
But  electrostatic  interaction  prevents  the  electron  distri- 
bution from  ever  being  completely  arbitrary,  and  it  is  from 
partial  spacial  coherence  that  so-called  collective  scat- 
tered radiation  is  observed.   The  degree  or  intensity  of 
coherence  is  dependent  on  the  ratio  Oi   between  the  wave- 
length and  the  characteristic  Debye  shielding  length  for 
the  plasma.   That  is, 

X v  (      KTS        "\^r 

d        =     zr      »      A^  =  {    *r   )2    ,    s    =   i,e 

4ff   XD   sin  |  D        V4ffne"   J 

(30) 

Where  X  is  the  wavelength  of  the  incident  light,  X  is  the 
Debye  length  and  6  the  scattering  angle.  For  a.  «  1  scat- 
tering is  incoherent;   for  &   »  1  spacial  coherence 
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determines  the  shape  of  the  scattered  spectrum.   Perhaps  a 
more  physical  and  qualitative  (if  not  mathematically  rigor- 
ous) understanding  may  be  had  of  this  concept  if  one  imag- 
ines the  Debye  spheres  as  integral  units.   Imagine  the  elec- 
trons as  tied  to  the  ions  they  shield.   Within  the  Debye 
sphere,  the  electron  distribution  is  arbitrary;   that  is, 
incident  of  wavelength  less  than  the  radius  of  the  Debye 
sphere  will  oscillate  electrons  within  that  sphere  in  dif- 
ferent phases,  whereas  incident  light  of  wavelength  greater 
than  the  Debye  sphere  will  oscillate  the  Debye  spheres,  and 
thus  all  electrons  within  them  in  unison,  in  coherence. 

The  incoherent  scattering  spectrum  is  a  direct  illus- 
tration of  the  electron  velocity  distribution,  as  will  now 
be  shown.   As  the  scattering  under  consideration  is  inco- 
herent, the  effect  of  a  single  electron  may  be  looked  at. 
It  has  been  seen  that  the  electric  field  of  light  scattered 
by  the  j    electron,  by  equations  (7)  and  (11)  is  propor- 
tional to  the  phase  factor  exp[i£0  t  +  ik-_r.].   Approximat- 
ing r_.  by  an  expression  valid  for  straight  line  electron 
paths, 

r_.(t)  =  r..(0)  +  v  t   ,  (31) 

the  scattered  field  has  the  frequency 

o>"  =  05  +  k-v .  (32) 

o : 

and  the  scattered  light  is  Doppler  shifted  by  k-v.  =  kv . 
(k  =  k  -k ' ) .   The  shift  is  two  fold;   that  is  the  electron 
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is  a  moving  source  to  the  scattered  light.   Thus  the  total 
shift  is  k*v  =  k"-v  -  k'-v.  F(v)  now,  is  the  normalized 
velocity  distribution  function  with  velocity  v  parallel  to 
k.   Corresponding  to  velocity  v,  the  frequency  is  Ui  =   kv; 
the  velocity  interval  dv  correspond  to  frequency  interval 
dco  =  kdv.   Thus,  the  number  of  electrons  with  velocity  v 
is  proportional  to  the  scattered  intensity  with  shift  kv. 
That  is, 


Fe(v,   -   ^g^  ,33, 


where  v  is  the  velocity  component  parallel  to  k,  I (kv,k)  is 
the  scattered  intensity  integrated  over  all  frequencies. 
Again,  this  is  only  valid  for  QL   «   1. 

The  topic  now  turns  to  the  more  complicated  theory  of 
longitudinal  plasma  waves.   The  attempt  here  is  to  calculate 
the  mean  thermal  value  of  the  square  of  the  space  time  de- 

pendent  Fourier  coefficients  |  n,  (£0)|  .   The  basic  reference 

~~  3 

for  this  calculation  is  the  paper  by  Salpeter.    The  ampli- 
tude squares  |  n,  (0J)|  "  will  be  solutions  of  the  collisionless 
Boltzmam equation,  i.e.  the  Vlassov  equation  (See  Appendix  A 
for  derivation  of  the  Boltzmann  equation) .   The  collision- 
less  Boltzmann  equation, 

-li  +  v-7  f  +  -  •  V    f  =  0  (34a) 

ot r    m   —v 


3 
E.  E.  Salpeter,  Phys.  Rev  12,  1528  (1960) 
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is  essentially  a  continuity  equation  for  representative 
points  in  electron  phase  space.   The  increase  in  time  in 
the  number  of  electron  points  in  phase  space  is  caused  by 
three  factors:   (1)  The  drift  in  points  in  position  coordi- 
nates due  to  inherent  velocity:   (2)   the  drift  in  points 
in  velocity  coordinates  due  to  electric  and  magnetic  forces, 
and;   (3)   the  spontaneous  velocity  change  without  position 
coordinate  due  to  collisions. 

In  the  purpose  of  calculating  |  n,  (u>)  |   the  spacial 

Fourier  coefficient  f   (t)  of  the  distribution  function 

qv 

f(r_,v,  t)  is  introduced.   (q  is  the  general  wave  vector). 
Corresponding  to  equation  (22)  , 

<-.(vdv) 
f  v(t)   =  I  expC-iq-r^t)]  ,  (34) 

~~~        J 

the  sum  being  taken  over  those  particles  with  velocities 
between  v+dv.   As  dv  -  0,  the  right  hand  side  has  one  term 
if  a  particle  of  velocity  v  exists  and  no  term  if  such  a 
particle  does  not  exist.   Thus  (34)  represents  a  fine  grain- 
ed distribution  function  so  described  because  the  density 
has  but  two  values :  zero  in  a  vacuum  and  infinite  at  the 
position  where  a  particle  is  represented  by  a  mathematical 
point.   Clearly,  then,  an  integral  over  the  velocity  dis- 
tribution becomes  identical  with  a  sum  over  all  the  parti- 
cles : 

f   (t)dv   =   >  exp[iq.r.  (t)]  (35) 

,00     _i J 

j=l 
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The  number  of  terms  in  (34)  is  equal  to  the  number  of  par- 
ticles with  velocities  in  v  and  v_  +  dv„   If  the  physical 
properties  that  are  influenced  only  slightly  by  the  micro- 
scopic details  of  the  distribution  function  are  calculated 
and  if  the  course  grained  (continuous)  distribution  is  in- 
dependent of  spacial  position,,  the  number  of  terms  in  the 
sum  is  the  most  important  aspect  of  (34) ,  and  may  be  approx- 
imated by 


f„„(t)  =  n(r,  t)F(v)  =  ^^~     £  n  Jt)exp(iq.r)      (36) 


F_ 
V     L>   ' 

q 


where  F(v)  is  the  normalized  coarse  grained  (continuous) 
velocity  distribution  and  n  (t)  is  the  Fourier  coefficient 
of  the  particle  density  given  by  (21)  and  (22).   As  the  no- 
tation in  (36)  implies ,  F ( v)  is  independent  of  time.   This 
of  course,  is  valid  only  if  Af/F  is  small  during  the  time 
1/Ao)  where  hue   is  the  practical  spectral  resolution. 

Without  regard  to  the  coarse  or  fine-grained  structure 
of  f  (t),  the  distribution  function  may  be  represented  for 
ions  and  electrons  as  a  Fourier  series s 

Fs(r,v, t)  =^     £f    (t)exp(iq.r),   s  =  i,e       (37) 

q   ~~ 

Similarly,    for    the    three    terms    in    the   Vlassov   equation, 
9f  .      r  *fsrM,(t) 


ts         1      V   atsc[v^ 

—  =  V     L  Bt  exp(iq.r) 


(38) 
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and  for  the  second  term 
3f 


s 
v 


=  ¥  '    Zsfsav(t)  exP(i£-£)  <39> 


—    Br     V    Lx   _±  sqv 

q 


In  the  third  term,  only  the  electrostatic  force  e  E_(r_,  t)  is 
considered  where  e   is  the  charge  of  the  particle  at  r*E_(r,t) 
is  written  as  a  Fourier  series, 


E(r,t)  =  ~  X^q(t)  exP(iq-r)  (40) 

For  -^—  equation  (36)  is  used,  where  F(v)  is  the  only  ve- 
locity dependent  quantity.   The  third  term  then  becomes, 

9        Sf    e    v  r  SF  (v) 


^[E  ,exp(iq'.r)2ns  „exp(iq".r) ] g- —   (41) 


m    ov       TT2  Lj  —a        F      -±  — '  L,   sq 
"isV   q 


s    —    m_V   q  —  q  — 


Combining  terms  in  the  double  sum  for  which 

q'  +  q"  =  q   ,   q"  =  q  -  q'  (42) 


^s   Bf     es    ^  „  SFs(v) 

— -  •  -c— ^  =     0  )  )  E  '  ,  n  ,    ,  N  exp(iq.r)  •  = (43) 

m    Bv       TT2  Lj   Z_,,— q'  s(q-q')  * x  -i  —'  Bv  v   ' 

s    -    mgV   q  q  ^    v^  % 


IV-      ,.    ,   SFs^ 


m  V 


2  Z,  ^Sg  exp(iq.r).  — ^ 


where 


^sq  "  esl\!  ns(q-q")  (44) 


23 


The  "Random  Phase"  approximation  is  now  introduced.   The 
order  of  magnitude  of  n   follows  from  (22)  .   Remember  that 


Lim   n   =  N  »  n   -  ,/W 


(45) 


n   is  the  sum  of  N  positive  and  negative  term  with  values 
arbitrarily  distributed  from  zero  to  one.   The  probability 
that  the  real  and  imaginary  parts  of  n   have  the  values 
£   and  T)  respectively  is 

q     q 


w 


(*a'V  =    2^N  eXP  [ 


q  q 


(Cg,+  y) 

2N 


(46) 


and  the  mean  amplitude  square  value  is 


00  CO 


<l  n   \> 

\  i      q   I    / 


d£   dT?    (4    2    +    T?  2)W(C    ,V   )    =    2N  (47) 

_»  si    q    q         q         q    q 


(46|  and  (47)  explain  the  importance  of  (45)  .   E   that  is, 
q  =  0  has  no  dominating  position  for  the  electric  field 
because  it  has  a  mean  value  in  a  neutral  plasma  of  zero. 
Thus  (44)  has  one  dominating  term  only. 
Writing  explicitly. 


r    =  e 
sq     s 


E  N  +  )  E  ,n 
L— a     La  —a   £ 


E  N  + 

— q     L— q'  s(q-qD) 

q'*q 


(48) 
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The  "random  phase"  approximation  assumes  that  the  sum  can 

4 
be  neglected  in  comparison  to  the  term  E  N.    Adopting  this 

approximation,  equation  (33)  finally  becomes 

1  V         rMsav  esns      SFS(Z)  -i 


q 


s    -± 


The  Fourier  representation  is  unambiguous,  and  thus  (49)  is 

valid  not  only  for  the  sum  of  all  general  wave  numbers  q 

but  also  for  every  single  q  .   Writing  (49)  explicitly  for 

electrons  and  multiplying  by  the  electron  charge  e,  the 

equation  for  the  Fourier  coefficient  <J  corresponding 

6 1  q*  v 

to   charge   density   a      : 

3  J      ev 


dcr  2 

eqy_  e   n  oF  (v) 

5T +  ivq   CT  +  E     .      s~=      =   0  (50) 

at  —  —     eqv_         m        — q         ov 


An   analytical    expression    for   E_  may  be   obtained   from 
Poisson's    equation 

7    .  E     =  4ff(J  =  4^(a  +cr.)  (51) 

—      —  e      i'  v 

Where  a  is  the  charge  density  created  by  electrons  and  ions 
As  previously  done,  let 

(J(r,t)  =  i£  <7  (t)exp(iq.r)  (52) 

q  - 


4 
Though  the  infinite  number  of  terms  in  the  sum  makes 
it  by  no  means  really  clear  that  the  approximation  is  al- 
lowable. 
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From    (52),     (51),    and    (40) 


V  iq.E      =   41T   Y  (T      ,    E      =    -4fficr     ^  (53) 

^     -  -£  ^     q        ~q  q  q2 

and  q  =  k  (as  in  equation  (12) ) 


ekv  ,   4ffn   ^     Bf  (v) 

+  i*v  cr  _   =  i  (  — S-  V  — |—  (54) 


St  ekv     v  mek  ^  k    Sv 

Where  v  and  E  are  the  components  of  v  and  E_  parallel  to  k 
(k  is  along  the  z  axis) .   The  Fourier  time  transform  of 
ff  ,   is  needed  in  the  calculation  of  the  spectral  distribu- 
tion;  this  is  written  as  in  equation  (20) 


qekv(t0)=J  aekv(t)expt-(iW+y)t]dt  (55) 

o 


Integrating  over  all  velocities, 


09 
P 


cx,  (o>)dv  =  QQV(a3)  =   cr  ,  (t)exp[-(iu>+y)t]dt      (56) 


,  ekv         ek 

00 


ek 
o 


q    (a?)  satisfies  partial  differential  equation  obtained  by- 
equating  the  transform  on  the  right  and  left  of  (54) .   The 
following  integral  is  needed 

-»   B<7 

— fj^  exp[-(iUH-y)  t]dt  =  -aekv(0)  +  (ia)+y)qekv(w)     (57) 
o 
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This  is  obtained  by  integrating  by  parts  and  using  (55) . 
Notice  that  the  contribution  from  the  upper  limit  a   ,   (°°) 
disappears  due  to  exp(-yt).   Using  (57)  in  (54)  gives 

.   47Tn  2.        SF  (v) 

q     w  (i^+y+ikv)-  crekv(o)+  i  Q  — — ;  o^m  -v      (58) 

ek  — 

or  multiplying  by  i      ~ 

r  4ffn   N       Bf  (v) 

iCTekv(°»-C^-)^(w'  ~V- 

ek  — 

q    (05)  =  ; (59) 

^ekvx  -05  -  kv  +  iy 


where 


Qk(0))  =  Qek(^)  +  Qik(w)  (60) 


The  charge  density  for  the  ions  is  obtained  similarly,  that 
is, 

,4ffn  2X       &F.  (v) 


i0..  (0)-f ^-Oa  i<4    — 5" 

lkv v  '    V  m  ,  y  ^k v  '    o^ 


ikv    V  m  ,   x  k      oy_ 
qikv(W)  =         -u>   -  kv  +iy  (61) 


f4*ne^     Tl  SFe^   1  *i.teh 


and  the  total  charge  density 

lek  '    JV   "*e   ~—   "1  "— 
qky_ ( W)  =qekv ( W)  +qikv ( ^  =  -05  -  kv  +i05 

(62) 
Integrating  (62)  over  all  velocities  gives  an  expression  for 
Q,(o?),  which  may  then  be  introduced  into  (59)  and  that 
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result  integrated  over  all  v.   The  final  result  would  be 
the  amplitude  necessary  for  scattering  Q  ,  ( 03)  =  -  en,  (a>). 
Since  the  integrals  may  not  be  actually  evaluated  without 
choosing  special  velocity  distributions,  the  following 
function  is  introduced 


v  p00  ^F  (v)   dv 
s  v     — 

s  =  i,e. 


_oo   ?v   -OJ-kv+iy  ' 

(63) 
Furthermore,  corresponding  to  (35)  the  velocity  integral  may 
be  calculated  as  a  sum  if  i(7,  (0)  is  chosen  as  a  fine-grained 
function. 
Then,  by  integration  of  (62) 


Qk(o>)  =  ie  I 


exp(-ikZ.)    exp(-ikz.) 


-to-kv.  . +iy   -0>-k  .  . +iy 


1  -  G  (o>)  -  G.  (0)) 


(64) 


Where  Z.,  z.,  v. .,  v   are  position  and  velocity  components 

y     y     id  ej     , 

of  ions  and  electrons  parallel  to  k  at  t  =  0. 

Equation  (64)  may  be  considered  the  basic  equation  of 
the  theory  of  longitudinal  plasma  waves.   From  this  equation 
follow  the  formulas  for  scattering,  dispersion  and  drift 
instabilities . 

The  next  step  is  to  employ  the  theory  of  density 

fluctuations  to  light  scattering  in  a  Plasma.   This  requires 

the  amplitude  squares  |  Qk(^)     and  |  Q  ,  (OJ)    ;   these 

represent  a  time  average  and  a  time  integral  over  all  micro- 
scopic states  assumed  by  the  plasma  between  t  =  0  and  t  =  — . 
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Substituting  (64)  into  (59)  and  (61)  and  integrating 
over  all  velocities  using  the  same  method  by  which  Qv(^)  was 
obtained,  |  Q  v^'   ma^  -be  calculated. 

-ie       f  (1-Gi(a?))^exP(-i^zi)  ^   £exp(-ikz,.)-| 
Qek(W)  =  1-G  (w)-G.  (a?)  L  -03-kv+iy   "^^  -uv-kv+iy  J 

(65) 

and 

n  ( ^  _      ie [  ( 1-G  ( (4  )  E  exp(-ikZ  )     E  exp f-ikz .)  -, 

Qek^  ;~  1-G  (05)-G  (05)  L     e     J J-_G.(05)3  X.1 

-     x   e^  ;  ^i^^  -to-kv+iy  ik  '   -OJ-kv+iyJ 

(66) 

The  amplitude  squares  |  Q  ,  (o>)|'   contain  double  sums  of  the 

form, 
E  exp ( ikz .       E  E  exp [ ik ( z  . -z . ) ] 

j  1 _  j  I    * 1 (c.n\ 

-w-kv  .+iy        (-ov-kv  .+iy)  (-os-kv  «+iy)  K      ' 

e  j  e  ~j  e  Jo 

as  well  as  double  products.   Ion  velocities  and  positions 

are  indexed  by  i,   whereas  electron  velocities  and  positions 

are  indexed  by  j.   Diagonal  terms  of  the  left  (the  single 

sum)  are  treated  the  same  as  off-diagonal  terms.   In  the 

double  sum,  however,  the  diagonal  terms  behave  differently; 

that  is,  the  space  dependent  phase  factor  vanishes  in  these 

terms  allowing  the  replacement  of  the  sum  with  an  integral 

over  the  velocity  distribution: 

N 


)  exp[ik(Z.-z.)  ]      p°°  F  (v)dv 


N      ,".w^../>!  (68) 


j  =1  (-0J-kv)^+  y^  '     J  _00(OH-kv)  ^+yz 


Since  the  denominator  of  the  right  hand  integral  depends 
only  on  v  ,  the  integral  is  over  the  one-dimensional 
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distribution  F  (v) ,  gotten  by  integrating  F  (v)  over  v   and 

v  .   Rigorous  evaluation  of  the  integral  is  at  best  tedious. 

It  may,  however,  be  approximated  more  readily.   The  main 

contribution  to  the  integral  comes  at  resonance  v  =  -  r-  . 

k 

At  resonance,  the  value  of  the  denominator  varies  rapidly 
while  the  numerator  varies  slowly,  permitting  it  to  come 
out  of  the  integral. 


lim 


F  (v)dv  p-       ,  ttf  (-r-) 

dv  e  k 


e x  '  „  ,  -as  . 

-  =  F  ( 


^  /  .  ,t  ^  2  ,   2     e   k  ' 
y->0   -  (w+kv)  +  y 


-°°  (o>fkv)2+  y2  yk 


(69) 


yielding  an  expression  for  the  diagonal  sum  dependent  on 
the  time  of  observation  —  .   This,  of  course,  is  what  is 

y 

expected,  since  |  Q  .  ( o)     determines  the  scattered  radia- 
tion during  T  =  —  . 

y 

The  off -diagonal  terms  behave  somewhat  differently.   It  will 

be  shown  that  this  sum  may  be  neglected  compared  with  the 

diagonal  terms.   The  sum  is  independent  of  v  and  may  be 

neglected  in  comparison  to  the  double  sum  (which  depends 

on  — )  for  small  y.      Consider  the  double  sum  over  a  factor 

expCik (z t-z . ) ]  and  assume  the  denominator  is  replaced  by  an 
*  3 

2 
appropriate  mean  value.   The  result  is  a  sum  over  N(N-1)-N 

positive  and  negative  terms.   From  the  same  "random  phase" 

the  sum  is  of  the  order  of  magnitude  N.   Hence, 


r-»  r-.      expL 

L  L    {-(x>- 


exp[ik (z  „-z.)  ] 


LZ± 


-    N 


00     00 

'   r  F    (v  . )  F    (vjdv.dv. 

e %    y     e      i'       i       I 

(-OJ-kv    .+iy)  (-o)-kv    ,-iy)        "  JJ^  (-ovkv  .+iy)  (-ui-kv.-iy) 
J  (70) 
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The  single  integrals  over  dv .  and  dv  have  linear  resonances, 

3  i 

which  by  symmetry  are  integrable;   the  numerator  remains  un- 
changed as  the  denominator  changes  sign.   Thus  the  resonance 
contribution  to  the  integral  is  zero.   Therefore  the  single 
integrals  are  independent  of  y   and  may  be  neglected  in  com- 
parison to  the  diagonal  terms  proportional  to  —  .   One  then 
obtains. 


<W*>'2    ,     ,..,2    »Nli-V^2Vx>+IV^|2Vir> 


eZ  K  7K        I  1-G    (0?)-G.  (0>)|  Z 


e  i 


Oik*""12        itn  I  l-Se(^)l2Fi(^)+  |G.(a,)|2Per-f) 


e 


yk      I  1-G    (W)-G.  (03)|  2 


(71) 


(72) 


QAUS)\2  W1VT      F    (^)     +    F.  (=£) 


e 


2  *      ll-G    (W)-G.(^|2  «"> 


Evaluation  of  these  amplitude  squares  of  course  depends  on 
calculating  G   and  G.  for  special  electron  and  ion  velocity 
distributions . 
Of  interest  is  the  shifted  Maxwellian  distribution; 

3 

Fs^>  -  (  «KtO  4  exP  ["  SET"  <2  +  V*]         <74> 

s  s 
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where  W   is  the  drift  velocity  of  particles  of  type  S, 
(S  =  i,e) 
Clearly, 


3F  (v) 
s  — 

5v~ 


-[ 


m 


^r-  (v  +  ws)J  fs(v) 


(75) 


and 


G  (O))  =  ka  ' 
s        s 


(v+W  )F  (v)dv 
x    s   s x 

(-w-kv+iy) 


(76) 


where  a.      is  given  by  (30)  and  G  (o>)  by  (63)  ,  and  where 

F  (v)  is  a  one-dimensional  velocity  distribution  in  v  ,  and 

W   is  the  z  component  of  W  . 

Introducing  the  Doppler  shifted  frequency,  one  has: 


a)   =  0J-  kW 

sw        s 


(77) 


Quantities  of  the  unshif ted  (no  drift  velocity)  distribu- 
tion will  now  be  superscripted  with  a  zero. 
Then, 


G  (oj)=  ka  ' 
s       s 


00  (v+W  )F°(v+W  )d(v+W  ) 

s   s  v    s'     x    s      o 

oo  [-o>  -k(v+w  )+iy]     :  Gs(V 
sw   %    s' 


(78) 


and 


F  (v)dv 
s 


2   2 
_oo  (oH-kv)  +y 


F  (v+W  )d(v+W  ) 
s v    s       s 

-00   [uj  +k(v+w  ]  +y 
sw   v    s 


ff   o 
— —  F 

yk  s 


-03 


SW 


(79) 


32 


Using    the   Doppler   shifted   frequencies 
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,  .2  2   n  2     2  ps  .  .  #««» 

"..■    =   3   k   Vs.s   =^2     '    s's      =   1'e  <80> 

s 

2 
where  v    ,  is  the  mean  thermal  velocity  square  of  the  par- 

Si  S 

tides  of  mass  m   and  temperature  T  ,  and  introducing  the 

s  s 

dimensionless  frequencies 


o>  o\ 

x    =  y    =  (81) 

05  J       0).  .  v   ' 

ee  li 


it   follows    that    (using    (74)), 

2      mh     r   |  1-Gi(y')|2  exp(-x'2)+(Ge(x')|    T?exp(-y'2)n 


»V«*-£t[: 


^ee    L  I  1-Ge(x')    -   Gi(y') 


(82) 


and, 


Y  03..  K2}        V..  '  05.  .T)        l2;  O)  ^J' 

ii  ii  ii  ee 


where 


"  =  Or;  Qsr ;  <84> 


Furthermore,    if    the    substitution 

u      =     W   +   v  (85) 
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is  made,    and   the    integrand  of    (79)    is    expanded  with    the 
factor    (-03-ku-iy)  ,    then, 


2     p°°       -"° 
G(x)    =  ko 


uF    (u)  (-03-ku-iy)  du  ,      . 

2       2  \&v) 

-00      (-OJ-ku)    +y 


The  imaginary  part  is  an  integral  of  the  form  of  (69) . 
Therefore, 

lim   I(x)  =  *%■   ax  exp(-x  )  (87) 

y-0 

Assuming  y  small  the  real  part  of  G(x)  may  be  expanded  in  a 
power  series,  which,  when  evaluated,  is  equal  to  zero  in  the 
limit.  Using  the  same  arguments  as  with  (70),  the  real  part 
goes  through  zero  at  resonance  and  changes  sign.  Again, 
integrated  over  symmetric  limits,  the  integral  goes  to  zero. 
From  F  (u) ,  let  <|  u|>  be-  the  mean  value  of  I  u| ;  expand- 
ing the  integrand  of  (86)  for  03  >  k  <C  I  u  |  >  in  the  form 

00 

-.  •         -0>-ku  1  .  -1     V/    i  \ v    Mv  /oo\ 

^  - — ;  ,2   2  =  rsasi =  "or  Z^-d    (-«)  <88> 

y-0      (-0!-ku)    +y  U=Q 

the  following  is  obtained : 

2 
lim   Re(x)  =  ^  (1  +  §  x  +  -^  x2  +  ---)  (89) 

y-0  ax 

For    03<k  <C|  u  |  y   ,    the    substitution, 


6  =  ku   -    w  ,      u  -  -  IfcSi-  .        du  =   -  f  (90) 
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is   made   and   F    (u)    is    expanded.      That   is, 


V   ,   m. 


P°(u>=F°(i±^)  .  P°(-  I)exp(.x2)  pzMl  (^L)  ,  ■> 

V=Q  k   KT 


(91) 


Using    this    to   obtain 


lim      Re(x)    =   -   a2[l-2x2(l   -  -|  x2   +jtx2+    •••)]  = 

y->0  u  J  Id  J 


a2[l-F(x)]  (92) 


where    f (x)    is    identified  with 

2       Tx  2 

f(x)    =   ttx   exp(-x    )         exp(t    )  dt  (93) 


because  this  function  has  the  same  series  expansion.  Thus, 
finally, 

Gs(xs')  =  -  as2[l-f(x^)-  iwh   xs'  exp(-x2,)]  s=i,e  x^  =  y' 

(94) 
What  has  been  presented  here  is  the  linear  theory  for 
scattering.   Appendix  B  gives  calculated  example  curves  of 
spectral  distributions  of  plasmas  with  drift  velocity  W 
(W  is  smaller  than  W,  ,  the  critical  velocity  at  which  non- 
linear drift  instabilities  are  excited)  W  is  in  multiples 

3  h 
of  the  electron  thermal  velocity  hr)   v   ;  at   is   01   ,    the 

2     ee         e 


35 


T 
parameter  identified  in  equation  (30)  and  e  =  — —  .   Each 

5  it  2tt 

spectrum  contains  11  curves   from  angles  0  =  0,   yr-  ,  ttt*  ••''tT 

between  W  and  k  with  W   =  |  W  |cos  0  .   The  curves  with  the  most 

it 
pronounced  side  maxima  correspond  to  0  =  ■=■  ,  the  only  symme- 

it 
trie  curve  corresponding  to  0  =  ■=■  with  the  horizontal  tan- 
gent at  y  =  0.   y  is  of  course  the  dimensionless  frequency 

— —  ,  and  extends  over  |  y|  ^3.5  representing  the  ion  line. 

ii 

Relative  maximum  intensity  I   is  given  for  each,  since  inten- 

m 

sities  vary  strongly.   Intensities  are  normalized  so  that 

|  1-G.  (y)|2exp(-x'2)  +  |  G  (x')|2  T)  exp(-y2) 

l(0J,k)  =  ; -5 (95) 

|  1-G  (x1)  -  Gi(y)  | 


5 
From  the  report  by  Theimer;   by  permission  of 
Professor  F.  Schwirzke. 
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III.   APPARATUS 

The  experimental  setup  itself  requires  a  myriad  of 
equipment.   The  schematic  of  figure  1  shows  the  equipment 
involved  and  its  layout.   The  steady-state  plasma  facility 
at  the  Naval  Postgraduate   School  is  ideal  for  a  systematic 
investigation  of  the  laser  light  scattering  plasma  diagnos- 
tic technique.   The  parameters  of  plasma  density  and  elec- 
tron temperature  may  be  varied  over  wide  range  in  the  sys- 
tem.  The  plasma,  whether  argon,  nitrogen,  or  helium,  is 
produced  through  a  hollow  tantalum  cathode  and  is  acceler- 
ated to  a  disc  anode  between  voltages  of  around  40  volts. 
The  beam  is  columnated  at  the  cathode  by  an  arrangement  of 
wound  electromagnets  producing  fields  of  several  hundred 
gauss.   The  plasma  density  may  be  controlled  by  varying  the 
gas  flow  at  a  needle  valve,  altering  the  background  pres- 
sure with  a  10  inch  diffusion  pump,  and  by  adjusting  the 
field  strengths  of  the  columnating  electromagnets.   Plasma, 
and  hence  electron,  densities^  as*  seen  in  the  theory,  have 
a  direct  bearing  on  the  scattered  light  intensities.   Pre- 
paration of  the  system  required  the  construction  and  place- 
ment of  perpendicular  window  ports  through  which  the  inci- 
dent light  entered  and  the  scattered  light  emerged.   The 
ports  were  approximately  a  half  inch  in  diameter  and  were 
mounted  on  removable  flanges,  as  shown  in  figure  6.   The 
plane  described  by  the  perpendicular  incident  and  scattered 
light  wave  vectors  was  arranged  in  such  a  way  as  to  intersect 
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the  plasma  jet  at  a  point  within  a  half  inch  of  the  tanta- 
lum cathode  outlet,  for  two  reasons.   One,  this  is  presum- 
ably where  the  beam  and  hence  the  electron  density  are  the 
greatest;   secondly,  little  is  known  about  the  plasma  para- 
meters at  the  cathode...  the  scattering  technique  should 
yield  much  information  concerning  these  parameters.   Further- 
more, laboratory  space  limitations  and  apparatus  configura- 
tion demanded  location  at  the  spot  chosen. 

The  incident  laser  radiation  was  provided  by  the  model 
130  high  peak  power  Q-switched  laser,  manufactured  by  Optics 
Technology,  Inc.   The  instrument  consisted  of  a  power 
supply/control  console  and  a  separate  optical  head,  which 
incorporated  a  Fabry-Perot  type  interferometer  and  laser 
cavity  mounted  on  a  granite  base.   The  Q-switch  is  a  device 
which  enables  the  laser  to  deliver  extremely  high  peak  powex 
pulses  varying  from  one  to  a  hundred  megawatts.   As  here  em- 
ployed, it  was  a  bleachable  liquid  system  contained  in  a 
small  cell  with  optically  flat  faces,  which  only  needed  to 
be  inserted  between  the  laser  ruby  and  the  interferometer 
for  operation.   Further  function  and  use  of  this  device 
will  be  discussed  later.   The  laser  crystal  is  a  production" 

quality,  60°  orientation,  three-eighth  by  three  inch  ruby 

o 
crystal,  emitting  radiation  at  6943  A.   The  crystal  is  op- 
tically pumped  by  two  linear  xenon  arc  lamps  with  three 
inch  arc  spacing  and  overall  diameter  of  eleven  mm,  mounted 
opposite  the  crystal  holders  in  a  polished  aluminum  cavity. 
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Cooling  is  accomplished  through  a  forced  draft,  water-cool- 
ed air  system.   The  rear  element  of  the  Fabry-Perot  inter- 
ferometer is  a  TIR  prism,  and  a  multi-element  glass  re- 
flector mounted  on  a  micrometer  driven  assembly  comprises 
the  front  element.   Alignment  of  the  interferometer  optics 
is  extremely  critical  for  lasing  action  to  occur.   Align- 
ment may  be  accomplished  in  the  following  way.   With  the 
top  of  the  laser  cavity  and  the  ruby  crystal  removed,  care- 
fully position  a  cw  laser  so  that  its  beam  is  in  line  with 
the  axis  of  the  ruby  crystal  holder  pins  and  passing  through 
the  center  of  the  exit  hole  as  shown. 


RUBY  LASER 


ruby   holders 
flash  lamps 


adjustable 

mount 


CW  LASER 


The  task  of  aligning  the  cw  laser  may  be  more  easily  accom- 
plished accurately  by  placing  a  pen  or  a  pencil  on  the  ruby 
pins,  holding  it  vertically,  and  adjusting  the  cw  laser 
until  the  beam  strikes  the  pencil  in  the  center  when  placed 
on  either  pin.   There  will  now  be  a  multitude  of  reflec- 
tions on  the  face  of  the  cw  laser,  one  from  the  TIR  prism 
and  one  from  each  of  the  glass  reflectors  in  the  microme- 
ter driven   interferometer.   Determine  which  one  of  these 
reflections  comes  from  the  ITR  prism  by  placing  a  sheet  of 
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paper  directly  in  front  of  the  prism;   Once  finding  the 
correct  reflection,  remove  the  paper  and  hand  adjust  the 
prism  until  its  reflection  is  as  close  to  the  cw  exit  aper- 
ture as  possible.   Then  tighten  its  screws  at  its  base  to 
secure  it.   Next  adjust  the  micrometer  screws  on  the  inter- 
ferometer until  the  multiple  reflections  made  by  its  ele- 
ments merge  on  or  near  the  cw  exit  aperture  and  prism  re- 
flection.  This  accomplished,  fringe  patterns  should  be 
observed,  and  the  ruby  laser  is  aligned. 

The  bleachable  dye  Q-switch  device  previously  mention- 
ed permits  the  formation  of  giant  pulses.   The  technique, 
usually  referred  to  as  passive  Q-switching  employs  a  small 
cell  containing  a  solution  of  one  of  several  metallo-organ- 
ic  compounds  known  as  phthalocyanines.   In  this  case,  a 
solution  of  cryptocyanine  in  methanol  is  used.   The  solu- 
tion strongly  absorbs  light  of  the  ruby  frequency,  and  this 

absorption  prevents  net  amplification  of  light  until  a 

3+ 
greater  proportion  of  Cr     as  have  been  pumped  into  the 

upper  energy  state  than  is  normal.   The  energy  input  con- 
tinues until  the  amplification  in  the  ruby  overcomes  losses 
due  to  absorption  in  the  cell.   To  this  point  the  laser 
smits  coherent  light  weakly.   A  small  amount  of  the  light 
bleaches  the  cryptocyanine  which  then  becomes  nearly  per- 
fectly transparent  to  the  ruby  light.   At  the  same  instant 
there  is  a  large  net  amplification  and  a  giant  pulse  con- 
taining all  the  stored  ruby  energy  develops  rapidly. 
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Following  the  pulse,  the  solution  returns  quickly  to  its 
absorbing  state.   Below  is  a  typical  phthalocyanine  absorp- 
tion curve : 


CO 


5000 


5500      6000 
wavelength 


6500    7000  7200   7500 


Bleaching  of  the  metal  cryptocyanines  is  effected  by  the 
saturation  of  the  absroption  at  the  ruby  frequency.   That 
is,  the  solution  is  in  the  bleached  condition  whrn  the  pop- 
ulation densities  in  the  upper  and  lower  states  of  the 
phthalocyanine  molecule  (s',S  in  the  diagram  below)  are 
roughly  equal.   Because  of  the  large  cross-sections  for 
capturing  ruby  light  only  small  numbers  of  cryptocyanine 
molecules  are  involved.   Consequently  saturation  occurs 
rapidly  and  relatively  few  quanta  from  the  laser  beam  are 
required. 
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excited  triplet   T' 
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singlet  state 
S' 


ground    state 
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triplet  state  T 


ENERGY  LEVEL 
DIAGRAM  FOR 
PHTHALOCYANINE 


The  solution  itself  is  prepared  from  the  concentrate 
by  diluting  one  part  cryptocyanine  to  ten  parts  methanol. 
The  cell  is  filled  almost  to  capacity  with  methanol  and 
nine  drops  of  the  1:10  solution  is  added  to  this.   This 
appears  to  give  optimum  Q-switching  without  requiring  too 
high  a  threshold. 

Calibration  of  the  laser  output,  that  is,  determining 
output  energies  and  pulse  widths  (and  hence  output  power) 
involves  the  use  of  a  PIN  photodiode  (See  Appendix  C)  and 
a  laser  radiometer.   The  HPA  4205  is  a  silicon  planar  PIN 

photodiode  with  a  response  time  of  less  than  one  nanosecond 

9 
and  a  frequency  response  from  dc  to  10   hz.   The  diode  was 

wired  according  to  the  following  diagram. 
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PHOTODiODE 


20  V 


The  photodiode  is  useful  in  conjunction  with  a  fast-rise 
oscilloscope,  such  as  a  tektronix  519,  for  recording  laser 
pulse  shapes. 

Laser  output  energy  was  measured  using  a  Westinghouse 
model  RN-1  Radiometer  which  is,  more  specifically  a  calor- 
imetric  bolometer.   It  is  designed  to  trap  the  energy  beam 
from  the  laser  in  a  bundle  of  fine,  insulated  copper  wire 
which  functions  both  as  a  calorimetric  m  ass  and  as  a  bolo- 
metric  element.   The  change  in  resistance  of  the  wire  ele- 
ment is  proportional  to  the  energy  absorbed  and  is  almost 
entirely  independent  of  the  distribution  of  the  energy  with- 
in the  unit.   The  bolometer  element  is  part  of  a  convention- 
al Wheats  tone  bridge  circuit  and  the  change  in  balance  of 
the  bridge  is  a  measure  of  beam  energy.   (See  figure  4.) 
When  used  in  conjunction  with  an  electronic  microvoltmeter, 
laser  energy  may  be  read  directly.   When  a  microvoltmeter 
reads  the  imbalance  in  microvolts,  the  reading  is  directly 
proportioned  to  the  incident  laser  pulse  energy, 

^observed  reading  x    f  ^-V  Laser  output  in  joules 

^  in  microvolts   ^ 

(94) 
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The  schematic  of  the  calibration  setup  is  as  follows: 


ruby  laser 


J- 


beamsplitter 


oscilloscope 

The  following  are  photographs  taken  on  a  Tektronix 
545A  oscilloscope  set  to  trigger  as  the  pulse  was  initiated, 
The  time  scale  on  the  three  non  Q-switched  photos  was  50^ 
s/div.,  and  on  the  fourth,  Q  switched  pulse  500  ns/div.  or 
0.5jlls/div.   This  gives  an  idea  of  the  effectiveness  of  in- 
creasing lasing  power  by  greatly  shortening  the  pulse 
width.   The  Q-switched  pulse  width  appears  to  be  of  the  or- 
der of  200ns.   The  average  pulse  length  of  the  non  Q-switch- 
ed pulse  appears  to  be  of  the  order  of  250fis.   Since  power 
is  energy  per  unit  time,  this  shows  an  increase  in  lasing 
power  by  a  factor  of  1250.   More  specifically,  a  typical 
low  power  pulse  will  put  out  approximately  0.5  joules.   For 
a  pulse  width  of  250(is  this  comes  out  to  two  kilowatts. 
The  Q-switched  pulse  width  of  0.2|isec.  the  power  is  now 
2.5  megawatts,   a  substantial  increase.   The  first  three 
photographs  show  3  pulses  each;   below  them  is  recorded 
the  average  energy  output  of  the  three  for  a  given  input 
energy  from  the  laser  power  pack.   This  was  the  same  for 
each  trace  on  a  given  photograph. 
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Non  Q-Switched  Pulses 
Sweep  Speed  =  50  jis/div, 


Input  1100  joules 
Output  0.19  joules 


Input  1225  joules 
Output  0.38  joules 
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Non-Q-Swi  tched 

Sweep   Speed   50    jLts/div. 


Input  1300  joules 
Output  0.55  joules 


For  these  first  3  photos,  a  plot  of  input  versus  output 
energy  is  made,  establishing  a  linear  relationship  above 
laser  threshold. 
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The  fourth  photo  is  included  as  a  fine  example  of  a 
Q-switched  laser  pulse.   Sweeptime  here  is  0.5  /is/division. 


The  scattered  light  is  analyzed  with  the  Baird  atomic 
modified  Eagle  Mount  three  meter  concave  Grating  Spectro- 
graph.  The  spectrograph  is  a  Rowland  circle  type  mounting 
for  which  the  focusing  conditions  are: 


r  =  RcosA 


A  =  Angle  of  incidence 
B  =  Angle  of  diffraction 
r   =   Slit  to  grating  distance 
r1  =  Camera  to  grating 

distance 
R  =  Radius  of  curvature  of 

the  concave  grating. 

N  =  Normal 

G  =  Grating 

S  =  Slit 

C  =  Camera 


47 


The  Rowland  circle  has  a  diameter  equal  to  the  radius 
of  curvature  of  the  grating  blank.  The  grating  is  mounted 
tangent  to  the  Rowland  circle. 

The  usual  grating  equation  applies.   That  is 

nX  =  d(sin  A  +  sin  B)  (95) 

n  =  Order  number 
>  =  Wavelength  in  Angstroms 
d  =  Grating  spacing  in  Angstroms 
A  and  B  are  of  opposite  sign  if  they  are  on  different 
sides  of  the  normal.   Spectrograph  controls  include  toggle 
switches  for  the  grating,  the  plateholder  (on  which  the 
photomultiplier  tubes  are  mounted)  and  focus. 

Rotating  the  grating  changes  the  angle  of  incidence 
A  (and  hence  B) .   The  Veeder-Root  counter  is  geared  so  as 
to  indicate  the  angle  of  incidence,  subtract  from  1000. 
The  focus  toggle  increases  or  decreases  the  plateholder  to 
grating  distance,  the  tilt  control  moves  the  plate  (mount- 
ing the  photomultipliers)  in  and  out  about  a  vertical  axis 
through  the  center  of  the  slit.   The  latter  two  controls 
are  not  as  important  when  using  the  photomultiplier  as  op- 
posed to  a  photographic  plate  when  a  careful  focus  is 
required. 

There  existed  no  calibration  chart  of  the  instrument 
to  predetermine  the  position  of  the  photomultipliers  for  a 
specific  wavelength.   This  must  be  accomplished  experiment- 
ally, using  a  calibration  source  of  known  wavelength,  for 
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example  a  helium  lamp.   Helium  has  two  lines  close  to  the 
6943  A  ruby  wavelength,  one  at  6730  A  and  another  at  7108  A. 
Sighting  visually  through  the  plate  aperture,  these  may  be 
readily  identifiedo   With  a  slit  cut  in  dark  cardboard,  the 
position  at  which  these  lines  strike  the  aperture  may  be 
marked.   Linearly  interpolating  between  the  two  wavelengths 
gives  a  rough  position  for  placement  of  the  photo  multi- 
plier.  A  fine  adjustment  may  be  made  when  the  whole  system 
is  aligned.   The  alignment  procedure  will  be  described 
shortly.   For  purposes  of  this  experiment  the  proper  grat- 
ing setting  was  found  to  be  991.2  (that  is,  with  the  angle 
of  incidence  A  equal  to  8.8  degrees.).   The  dispersion  may 
also  be  calculated  with,  some  ease.   This  will  become  of 
great  use  when  actual  spectral  distributions  are  analyzed 
carefully.   This  relation  may  be  derived  from  the  basic 
grating  equation. 

n\   =  d(sin  A  +  sin  B)  (95) 

\   =   —(sin  A  +  sin  E) 
n 

letting   d  be    a   constant 

(||)  £cos    B  (96) 

d 

or  with  A   constant 


dX  d 

_.  cos    B 

dB  n 
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=   rdB 


(97) 

where  r  is  the  focal  distance  and 
is  the  increment  along  circumference 
for  dB  increment. 


dB 

1 

di 

r 

dX 

dX 

dB 

&l 

dB 

d£ 

(98) 


-77   =   —  cos  B  •  — 
dJfc     n         r 


(99) 


but 


r   =   R  cos  B 


(100) 


hence 


dX 


d   cos  B 
nR  cos  B 


D  =  actual  grating  spacing, 
d  =  effective  grating  spacing. 


dX 
di 


_d^ 
nR 


(101) 


d   =   D  cos  B 


(102) 


therefore 


dX      D  cos  B 
d£  nR 


(103) 
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Numerical  calculations  for  dispersion  may  now  be  made,  let- 
ting D  be  the  spacing  for  15000  groove/inch  grating  and 

o 
R  =  3  meters  at  the  6943A  wavelength  and  an  angle  of  inci- 
dence A  (determined  experimentally)  of  8.8  . 


1 .   ,      2  54mm    1Q8A     -ICn-5o£ 

D  =  t  c^o^   inches  x  — r- — r-  x  — =  16933A 

15000  inch     cm 


for  the  first  order  spectra 


X  =  d(sin  A  +  sin  B)        A  =  8.8° 


*  •      TV 

sin  B  =  —  -  sin  A 
d 


6943      .   Q  Oo 
16933  "  Sin  8"8 


41003  -  .15299 


25704 


cos  B   =  .96638 


d\    _   D  cos  B  ..  (16933)  (.96638) 
dl  R  3  x  103 


16369        _  A-    °, 
— — -r      =   5.45  A/mm 

3  x  10 
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The  photomultiplier  tubes  were  RCA  7102  tubes  with 
sensitivity  and  response  curves  given  below. 
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supply  volts 

The  tubes  were  operated  from  a  number  of  power  sup- 
plies, the  best  of  which  appeared  to  be  the  Power  Designs 
Model  HV  1565  regulated  dc  power  supply.  Optimum  signal 
to  boise  ratio  was  obtained  at  a  voltage  setting  between 
1450-1490  volts  (actual  voltage  as  read  on  the  voltmeter 
was  closer  to  1350  volts) .   The  output  signal  was  read 
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across  about  a  five  megohm  impedance  for  greater  signal 
strength.   This  was  of  course  at  some  sacrifice  to  response 
time,  as  will  be  seen  later. 

Before  discussing  the  important  alignment  procedure, 
the  optics  used  to  transfer  the  scattered  light  from  the 
plasma  to  the  spectrograph  will  be  described.   The  arrange- 
ment was  a  periscope  type  as  shown  in  figures  1  and  5.   Two 
converging  lenses  were  employed,   one  focusing  at  the  plasma, 
the  other  refocusing  the  plasma  light  at  the  exit  of  the 
periscope.   A  quartz  condensing  lens  on  the  optical  bench 
of  the  spectrograph  refocused  the  light  on  the  slit.   An 
arrangement  of  rigid  cardboard  packing  tubes  and  aluminum 
mirror  holders  formed  the  remainder  of  the  device.   The  po- 
sition of  the  periscope  exit  (on  the  plane  of  which  the 
plasma  beam  is  focused)  is  determined  by  using  the  grating 
lamps  and  quartz  condensing  lens  in  the  following  way.   The 
grating  lamps  are  hung  over  the  front  face  of  the  grating 
holder  and  switched  on  from  a  control  on  the  external  power 
supply.   The  quartz  condensing  lens  is  placed  in  its  holder 
on  the  optical  bench  and  left  free  to  move  along  the  length 
of  the  bench.   The  image  of  the  grating  lamps  may  be  observ- 
ed by  holding  a  white  sheet  of  paper  behind  the  lens.   Mov- 
ing the  paper  away  from  the  lens  will  eventually  focus  the 
grating  lamps  in  two  distinct  images.   Moving  the  paper 
further  away,  the  images  will  merge,  and  it  is  at  that 
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point  that  the  source,  in  this  case  the  exit  of  the  peri- 
scope, is  to  be  placed.  The  lense  is  now  secured  and  the 
periscope  is  mounted  in  its  proper  position. 

Because  of  the  distance  it  was  necessary  to  convey 
the  lights  and  because  of  the  small  image  of  the  lenses, 
alignment  was  extremely  critical  and  somewhat  difficult  to 
obtain.   Accurate  laser-plasma-periscope-spec trograph  align- 
ment may  be  obtained  while  the  plasma  system  is  at  atmos- 
pheric pressure  and  the  tantalum  cathodes  are  being  replaced. 
With  the  cathode  flange  removed,  remove  also  the  left  win- 
dow flange.   Place  a  wooden  slate  across  two  stationary 
bolts  of  the  cathode  flange  in  such  a  way  as  to  have  the 
slat  horizontal  and  just  below  the  center  of  the  vacuum 
chamber  (and  hence  below  the  plasma  beam  axis) .   The  care- 
fully insert  a  rigid  aluminum  or  stainless  steel  rod  down 
the  chamber,  supporting  one  end  inside  the  annular  disk 
anode  and  the  other  end  on  the  wooden  slat.   The  slat 
should  be  of  such  thickness  that  the  axis  of  the  rod  it 
supports  is  parallel  to  and  just  below  the  chamber  axis. 
This  arrangement  will  now  support  the  mirror  proble  to  be 
inserted  from  the  left  window  opening;   the  probe  will  act 
as  an  artificial  "scatterer"  so  that  proper  alignment  both 
of  the  optical  transfer  system  and  of  the  grating  may  be  ob- 
tained.  The  mirror  probe  consists  of  a  straight  wooden 

o  o 

dowel  with  one  end  cut  to  45  .   To  the  45   cut  is  applied 
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a  small  cutout  piece  of  polished  stainless  steel  to  act  as 
a  reflector.   This  may  be  glued  on  and  a  mark  may  then  be 
made  on  the  opposite  end  of  the  dowel  to  indicate  the  pro- 
per vertical  mirror  position.   The  dowel  is  then  inserted, 
one  end  supported  on  the  rod  arrangement  described  above 
and  the  other  end  on  a  notched  cardboard  edge,  taped  to 
the  window  opening.   The  supporting  rod  may  vary  the  position 
of  the  reflector  along  the  chamber  axis  without  significantly 
changing  its  position  perpendicular  to  the  axis;   pushing 
or  pulling  the  mirror  rod  itself,  varies  its  position  per- 
pendicular to  the  chamber  axis. 

The  next  step  is  to  place  a  strong  light  source  such 
as  a  mercury  lamp  before . the  grating  (inside  the  specto- 
graph)  and  aimed  back  toward  the  slit.   The  Mercury  beam  may- 
be focused  at  the  slit  using  its  own  lens  for  more  intensi- 
ty.  With  the  laboratory  darkened,  carefully  adjust  the  per- 
iscope arrangement  until  the  light  from  the  Mercury  falls 
on  the  supporting  rod,  which  has  been  located  parallel  to 
the  vacuum  chamber  center  axis.   This  done,  the  mirror  prop 
may  be  adjusted  so  that  the  light  falls  on  its  reflector. 
From  the  end  of  the  mirror  proble,  and  peering  into  the 
chamber  through  the  left  window  mount  the  probe  may  be  ad- 
justed carefully  so  that  its  light  falls  on  the  half  inch 
laser  window.   The  axis  of  the  mirror  should  be  vertical, 
that  is,  reflecting  straight  up.   With  the  Mercury  light, 
(which  may  be  easily  discerned)    shining  on  the  half  inch 
laser  entrance  window,  the  system,  save  the  laser,  should 
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now  be  in  perfect  alignment.   The  next  problem  is  to  align 
the  laser  so  that  its  artificially  "scattered"  light  may  be 
seen  by  the  grating;  and  when  this  is  done  the  grating  must 
be  rotated  until  the  laser  frequency  is  centered  on  photo- 
multiplier  tubes  at  the  spectrograph  exit. 

The  laser  is  mounted  on   an  aluminum  vertically  adjust- 
able stand,  the  laser  itself  containing  tilt  adjustment  on 
its  base.   Alignment  along  the  horizontal  plane  may  be  made 
more  easily  with  a  cw  laser  by  raising  the  axis  of  the  cw 
laser,  the  pulse  laser  may  replace  along  the  same  axis  and 
trial  and  error  adjustment  made  until  alignment  in  the  ver- 
tical plane  is  achieved  at  will.   Once  aligned  the  laser 
pulse  should  be  seen  in  the  slit  when  fired. 

The  grating  may  then  be  rotated  to  the  proper  wave- 
length.  The  phototubes  and  grating  have  already  been  rough- 
ly adjusted  according  to  the  helium  line  procedure  described 
above.   First  the  PIN  photodiode  is  positioned  with  an  arran- 
gement of  beamsplitters  (microscope  slides)  to  analyze  the 
laser  pulse.   The  signal  from  the  PIN  photodiode  is  then  fed 
to  the  upper  trace  of  a  dual  beam  oscilloscope,  and  the  sig- 
nal from  the  photomultiplier  to  the  lower  trace.   Upper 
trace  gain  is  lowered  because  of  the  intensity  of  the  laser 
pulse,  but  the  oscilloscope  arrangement  is  as  shown  in 
figure  14.   The  grating  is  rotated  by  0.1  degree  at  a  time 
in  order  to  find  the  optimum  position.   Below  are  the  traces 
illustrating  how  this  was  accomplished.   The  upper  trace  is 
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the  pm  tube  signal  and  the  lower  trace  is  the  laser  pulse. 
Beside  each  is  the  respective  grating  setting  and  its  cor- 
responding angle  of  incidence. 


o, 


991.2  (8.8") 


991.1  (8.9°) 


991.0  (9.0°) 


58 


991.5    (8.8U) 


991.3     (8.7°) 


991.2    (8.8°) 


Sweepspeed  was  50  /is/div.   Both  traces  were  triggered  by 
the  laser  pulse,  the  lower  trace  of  each  group.   The  laser 
trace  gain  was  0.5  V/div.  and  the  pm  tube  trace  25  V/div. 
All  were  non  Q-switched  pulses.   Note  the  difference  in 
pulse  shapes.   As  was  mentioned  before,  this  was  due  to 
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the  large  resistor  in  the  signal  line  (about  five  megohms) 
which  considerably  lengthened  the  rise  time  of  the  tubes. 
This  of  course,  is  no  consequence  since  relative  intensi- 
ties are  what  is  important  here  and  in  analyzing  the  spec- 
tral distribution  of  the  scattered  light. 

The  system,  now  constructed  and  aligned,  may  be  em- 
ployed in  the  attempt  to  discern  a  scattered  light  signal 
and  analyze  the  spectral  distribution  of  the  scattered 
light.   The  few  attempts  thus  far  have  been  unsuccessful. 
Some  combinations  of  detecting  tubes  and  oscilloscope  am- 
plifiers have  been  tried  to  increase  the  signal  to  noise 
ratio.   Photographs  have  been  taken,  an  example  of  which 
is  below. 
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The  photograph  on  the  bottom  was  triggered  at  the 
scope  with  no  laser  pulse  in  order  to  provide  a  noise  refer- 
ence from  the  pm  tube.   Tube  voltage  was  set  at  1490  (Ac- 
tually 1390  on  the  voltmeter)  .   Sweep  speed  was  50  jls/div. 
and  the  background  pressure  in  the  vacuum  chamber  was 
8  x  10    mm  Hg.   Oscilloscope  gain  here  was  0.05  volts/div. , 
as  shown  in  figure  14.   If  a  signal  is  present,  as  a  few  of 
the  traces  such  as  the  one  above  seem  to  indicate,  it  is 
virutally  impossible  to  distinguish  it  from  the  noise. 
Possible  improvements  will  be  discussed  in  the  conclusion. 

Attempts  at  focusing  the  laser  beam  seem  to  have  no 
effect,  but  this  is  almost  certainly  due  to  the  poor  quali- 
ty of  lenses  available. 
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IV.   CONCLUSIONS 

The  scattering  apparatus  has  been  erected  and  aligned. 
The  attempts  which  have  thus  far  been  made,  to  discern  the 
scattered  light  and  analyze  its  spectrum  have  been  unsuccess- 
ful.  The  plasma  parameters  specifically  the  density,  affect 
scattering  intensities  as  much  as  many  varied  instrument 
parameters  affect  the  experimenter's  ability  to  detect  the 
light.   All  oscilloscopes,  save  a  very  few  are  in  dire  need 
of  realignment  and  recalibration,  as  are  the  plug-in  units 
that  accompany  them.   When  the  scattered  light  is  finally 
detected  up  to  seven  photomultipliers  will  be  needed  to  si- 
multaneously analyze  the  spectral  distribution  of  a  single 
pulse;   knowing  the  dispersion  of  the  grating,  the  photo- 
multipliers,  connected  to  the  spectrograph  plate  with  short 
lengths  pipe,  may  be  accurately  spaced  around  the  central 
frequency  to  observe  the  whole  spectral  distribution  of  the 
scattered  pulse.   These  should  be  of  the  RCA  7102  type. 
The  7102  tubes  now  available  appear  to  have  lower  signal 
to  noise  ratios  at  sensitivities  required  that  one  would 
expect.   Improvement  of  that  ratio  could  be  of  significant 
importance  in  detecting  the  scattered  light.   Increasing  the 
plasma  density  by  purging  more  gas  into  the  plasma  system, 
increasing  the  field  on  the  cathode  and  anode  magnets  and 
by  controlling  the  main  diffusion  pump  may  well  be  the 
answer  sought.   Furthermore,  the  laser  system  available, 
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while  powerful,  may  not  be  quite  what  is  needed,  though 
focusing  of  its  beam  through  a  lens  of  better  quality  may 
help.   A  more  energetic  system  is  now  on  order.   It  is 
interesting  to  note  previous  experimental  attempts  to  study 
at  light  scattering  and  the  problems  mentioned  in  these 
efforts.   One  group,  Thompson  and  Fiocco  at  MIT  measure  the 


light  of  a  laser  scattered  from  an  electron  beam.   The  elec- 

-18, 


9    -3 

tron  density  was  about  5  x  10   cm   ,  and  the  ratio  of 


scattered  to  incident  intensity  was  of  the  order  of  10 

Applying  their  technique  to  a  low  electron  density  plasma, 

13    -3 
n  w  10    cm   ,  they  improved  the  ratio  of  scattered  to  in- 

-15 
cident  power  to  10 

The  best  scattered  to  incident  ratio  came  from  E. 

Funfer,  B.  Kronast,  and  H.  J.  Kunze,  who  measured  light 

17 
scattered  from  a  high  density  pinch  plasma  (n  «=s  10   )  . 

-11 
Their  ratio  came  to  10 

Distinguishing  the  scattered  light  from  background 
plasma  radiation  is  a  considerable  problem.   One  possible 
solution  that  may  be  tried  is  to  eliminate  the  plasma  ra- 
diation with  polarizers,  since  at  90   the  scattered  light 
will  be  polarized,  whereas  the  plasma  background  radiation 
will  not.   One  other  wasy  of  eliminating  plasma  radiation 
would  be  with  the  use  of  a  differential  amplifier  which  will 
differentiate  the  scattered  signal  from  a  pair  of  pm  tubes 
one  seeing  the  scattered  light  and  radiation,  the  other  only 
radiation.   This  technique  was  attempted,  but  the  character- 
istics of  the  two  pm  tubes  available  were  not  close  enough 
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for  comparison.   Furthermore,  the  differential  amplifier 
plug-in  units  on  the  oscilloscopes  were  not  functioning 
properly. 

The  above  has  summarized  experimental  progress  to 
date.   This  progress  should  enable  future  experimenters  to 
resume  development  of  the  experimental  technique  of  laser 
light  scattering  as  a  plasma  diagnostic. 
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APPENDIX  A.   DERIVATION  OF  THE  BOLTZMANN  EQUATION6 


Each  electron  is  represented  by  a  point  in  six  dimen- 
sional electron  phase  space.  The  coordinates  of  the  phase 
space  are  x,  y,  ,  zr  ,u  .,  •  u.,  u  .   Because  each  electron  has 

r-  -L  j£»      y"      2 

a  velocity,  the  representative  point  for  each  electron  must 
change  its  position  coordinates  with  time.   In  addition, 
any  change  in  the  velocity  of  the  electron  will  also  cause 
the  representative  point  to  move  in  electron  phase  space. 
The  velocity  of  an  electron  will  change  if  a  force  acts  on 
the  electron  or  if  the  electron  collides  with  any  other 
particle. 

Consider  a  fixed  differential  volume  in  electron  phase 
space  dV  =  dx  dy  dz  du  du   du   (the  "electron"  subscript 
will  be  omitted) .   An  expression  will  be  written  for  the 
number  of  representative  points  in  this  volume  as  a  func- 
tion of  r_,  u,  and  _t;   that  is,  it  will  be  a  "continuity" 
equation  for  the  representative  points  in  electron  phase 

space.   The  increase  with  time  in  the  number  of  electron 

3(f  dV) 
"points"  in  dV,  — yr is  caused  by  three  factors: 


1.  The  drift  of  points  in  position  coordinates  due 
to  inherent  velocity. 

2.  The  drift  of  points  in  velocity  coordinates  due 
to  electric  and  magnetic  forces. 


6 
From  Uman,  p. 34-36 
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3.   The  spontaneous  velocity  change  without  position 
coordinate  change  due  to  collisions. 

When  these  three  terms  are  evaluated,  the  Boltzmann  equa- 
tion for  electrons  may  be  written.  Note  that  the  effects 
of  external  gravitational  forces  and  inteparticle  nuclear 
forces  are  omitted. 

First  determined  is  the  increase  with  time  of  the 
number  of  representative  points  in  dV  due  to  drift  in  po- 
sition coordinate  space.   The  number  of  representative 
points  which  cross  the  surfaces  of  dV  in  time  de  because 
of  inherent  velocity  must  be  calculated.   The  velocity  of 
a  point  in  the  x  direction  at  x,  y,  z,  u,  u  ,  u   is  u  . 
The  number  of  representative  points  in  a  volume  of  "base 

area"  dy  dz  du   du   du   and  height  u  dt  will  enter  dV  in 
x   y   z        3    x 

time  dc  through  a  face  normal  to  the  x  axis  at  x,  y,  z, 

u  ,  u  ,  u  .   This  number  is 
x    y    z 


T  a(cTN)  "I  dfc  =  f  (r  u  t)u  dtdydzdu  du  du  (1) 

L    dt   Jv      e  —  —    x   J     x  y  z 


The  change  in  the  number  of  representative  points  in  dV  in 
time  dt  due  to  the  exit  of  points  through  the  equivalent 
surface  at  x+dx  is  given  by 


r  ii|!sL]      dt  = .  rf  u  +  3(f~u-> 

L         at       -J  y+dx  Lex 


J  x+dx  Lex  dx 


e   x' 

dx 


dt   dydzdu   du    du 
J  x     y      z 


(2) 
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Therefore  the  net  rate  of  change  of  points  in  dV  due  to 
velocity  in  the  x  direction  is 


3(d6N) 
It 


x 


J  x+dx 


3(u  f  ) 

x  e 

Sx 


dV 


(3) 


Thus  the  total  rate  of  change  of  points  in  dV  due  to  in- 
herent velocity  may  be  written  as 


[ 


3(d6N) 

at 


-  (u   -  Vrf  ) dV 
v— e       e' 


(4) 


Now  consider  the  drift  of  representative  points  in  velocity 

coordinates  due  to  electric  and  magnetic  forces.   The 

"velocity"  of  a  representative  point  in  the  u   direction  in 

+   dux  N     +   Fex  . 
velocity  coordinates  is  \   ~^~~  J    ~  \      m       )  •   Using  this 


dt 


m, 


"velocity"  the  net  influx  of  points  through  the  two  surfaces 
perpendicular  to  u  may  be  calculated: 


r  a(d6N)  |         r  B(d6N)^ 

L      st     Jux  +  L      bt    J, 


ux+du 


L  Bu   v  m     e  J . 
x     e 


dV 


(5) 

The  extension  of  this  expression  to  include  the  flow  through 
the  surfaces  perpendicular  to  u   and  u   follow 


B(d6N) 

at 


i-i 


F 

— € 

m 


Vu  f. 


dV 


(6) 
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At  this  point,  the  time  rate  of  increase  of  represen- 
tative points  in  dV  is  equated  to  the  influx  from  the  two 
drifts  (neglecting  collisions)  and  the  Boltzmann  equation 
is  obtained: 


C -si ) dv  -  -  (s.  •  VrO dv  -  C  ?  -™Odv 

e 

or  (7) 


e 


+  u      •    7rf      +  —    .    7uf      =   0 


?>t       —  e  em  e 

e 
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APPENDIX   B 
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APPENDIX    C 


4mc&$ 


ULTRAFAST    LOW    NOISE 
SILICON    PIN    PHOTODIODE 

The  hpa  4205  is  a  silicon  planar  PIN  photodiode  for 
etection  of  light  in  the  visible  and  near  infrared  regions, 
esponse  to  blue  and  violet  light  is  unusually  good  for  a 
ery  low,  dark  current  silicon  photodiode. 

Speed  of  response  of  this  diode  is  less  than  1  nanosecond, 
aser  pulses  shorter  than  0.1  nanosecond  may  be  observed, 
he  frequency  response  extends  from  DC  to  1  gc. 

The  low  dark  current  of  less  than  150  picoamperes  en- 
)les  detection  of  very  low  light  levels.  The  quantum  detec- 
on  efficiency  is  constant  over  six  decades  of  light  intensity, 
roving  excellent  dynamic  range. 

An  extremely  fast  photon  coupled  pair  can  be  obtained 
!  using  the  hpa  4205  in  conjunction  with  the  hpa  4107 
nfrared  light  source  in  an  identical  package),  for  con- 
:nient  use  in  card  and  tape  readers,  encoders,  and  similar 
jplications. 
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OPTICAL  CHARACTERISTICS 

Test  Conditions 

Response  at  7700  A 0.75  electrons/photon V  =  —  20  v 

0.25  MA/mW/cm2 RL  =  <  1  Mn 

0.125  /iA/pW 

Sensitive  Area 0.5  x  10~3  cm2 

0.010"  diameter 

Speed  of  Response Less  than  1  nanosecond 


ELECTRICAL  CHARACTERISTICS  (25°C) 

Test  Conditions 

Dark  Current   1 50  picoamps  maximum VR  =  —  10  volts 

Forward  Current ;...  50  mA  minimum VF  =  1  volt 

Capacitance    0.7  pf  typical VR  =  —  10  V 


MAXIMUM  RATINGS 

Steady  Reverse  Voltage 50  volts 

Peak  Inverse  Voltage  200  volts  (see  Note) 

Power  Dissipation 50  mW     T  ==  25°C 


Note:  Exceeding  the  PIV  will  cause  permanent  damage  to 
the  diode.  Forward  current  is  harmless  to  the  diode,  within 
the  power  dissipation  limit.  For  optimum  performance,  the 
diode  should  be  reverse  biased  at  between  5  and  20  volts. 
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Figure  1.   Spectral  Response 


TYPICAL     DARK    CURRENT  (NANO-AMPS) 
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Figure  3.   Typical  Capacitance  Variation  with  Reverse  Bias 
(f  =  1  mc) 
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Figure  2.   Dependence  of  Dark  Current  on  Temperature 


Figure  4.   Relative  Directional  Sensitivity 
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Figure  5.  Schematic  Diagram. 
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Figure  6.  Relative  Directional  Sensitivity  of  the  PIN  Photo- 
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